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ABSTRACT 

In  this  pax>«r  we  analyze  the  steady-state  bifurcations  from  the  trivial 
solution  of  the  reaction-diffusion  equations  associated  to  a model  chemical 
reaction,  the  so-called  Brusselator.  The  present  analysis  concentrates 
on  the  case  when  the  first  bifurcation  is  from  a double  eigenvalue.  The 
depiendence  of  the  bifurcation  diagrams  on  various  parameters  and  perturba- 
tions is  analyzed.  The  results  of  Technical  Suttmary  Report  #1844  are 
invoked  to  show  that  further  complications  in  the  model  would  not  lead  to 
new  behavior. 
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SIGNIFTCANCE  AND  EXPLANATION 


It  is  i;»' 1 ir.*v«d  that  [■atrtcrn  formation  by  bioloqi"a]  — <'  bf,- 

related  to  re^fatod  bifurcations  of  the  reaction-diffusion  '-n'  govern- 

in'! the  concentrations  of  the  various  chemicals  present  in  *i.i'  tem.  The 
full  com£)lexity  that  these  equations  can  exhibit  is  far  fr^m  understood. 

Here  we  take  an  initial  step  in  this  direction  by  analyzin':  th'--  tv  *'ircTt  Lons 
of  a model  chemical  reaction,  the  so-called  tri-moleoui  ar  m'rdr  i .r'  1,'^fever 
and  Prigogine,  at  and  near  a double  eigenvalue.  The  dopen  hTioc  this 
system  on  various  parameters  and  its  sensitivity  to  per tur!  u t iorr.  is  also 
stud ied . 

This  paper  differs  from  earlier  work  in  several  points.  First  we  apply 
the  notion  of  universal  unfoldings  (developed  in  MRC  Technir-al  E in'-  iry  Heiiort 
#1844)  to  show  that  no  qualitatively  new  behavior  would  result  frur  tnu  addi- 
tion of  further  perturbations  in  the  model.  Second,  we  sh'iw  that  r ymmr-try 
plays  a dominant  role  in  determining  the  behavior  of  the  system,  a fact  that 
seems  to  have  escaped  notice.  Finally  we  have  been  able  to  attaiii  a higher 
level  of  rigor  than  previous  work  on  double  eigenvalues. 
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The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  refiort. 
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SI-  Discusslori  of  thn  i-robljnm 

In  this  pdi-ir  wf  consider  the  onc-d  imens  iond  1 rear  t ion-d  i f f ns  inn  rquations  associated 
to  the  so  caUnd  tr i -molncu lar  model  of  hefever  and  I'riqoqine  if.|,  less  formally  known 
as  the  "Dtussolator" . For  the  i>ardmeter  ranqe  in  which  we  are  interested,  this  system 
exhibits  bifurcation  from  a spatially  and  tem[«raHy  homoqeneous  solution  into  steady- 
state,  Sfiatially  inhomoqeneous  solutions.  The  first  bifurcation  may  be  from  either 
a simple  or  a double  eigenvalue-  A number  of  authors  11,11  have  discussed  the  case 
of  a simple  eigenvalue,  but  rather  less  seems  to  be  known  about  the  non-simple  case. 

Mere  we  analyze  the  bifurcations  of  this  system  at  the  double  eigenvalues  as  an  application 
of  the  theory  developed  in  12,31-  We  obtain  a rather  comfilete  classification  of  the 
[jossible  bifurcation  diagrams  in  the  vicinity  of  such  {Xiints.  Some  of  our  results 
were  obtained  earlier  by  Keener  15a I by  less  rigorous  methods  - his  results  are  compared 
with  ours  at  the  end  of  51. 

The  relevant  equations  for  this  model  are 

^ = D + X^Y  - (B  + IIX  t A 

^ 3f.2 

(1.11 

av  3^Y  2 

_ = X Y t BX 

subject  to  boundary  conditions  of  Dirichlet  type 

(1-2)  X(0)  = X(TI)  =*  A,  y(0)  ^ Y(tt)  « B/A  . 

Here  the  unknown  functions  X and  Y are  chemical  concentrations,  A and  B are 
constant,  externally  controlled  concentrations,  and  are  diffusion  coefficients. 

B plays  the  role  of  the  bifurcation  parameter;  that  is,  we  are  interested  in  the 
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bifurcation  of  ntw  solutions  of  (l.M,  (1.2)  from  the  trivial  solution 
( 1 . i ) X =■  A , 'f  = B/A 

us  B IS  increased.  To  facilitate  the  analysis  we  define  incremental  variables 

u = X-A,  v = y-  B/A  , 

for  which  we  use  the  vector  notation  w = (u,v).  These  variables  satisfy  the  equation 


(1.4) 
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B<jth  time  independent  and  time  periodic  .solution.s  of  (1.4)  can  bifurcate  from  the 

zero  solution  of  this  equation,  dependinq  on  the  various  parameters.  However  if 

is  rather  larger  than  say  ^2/0^  is  at  least  3,  then  the  first  bifurcation 

gives  a time  independent  solution.  (.See  (11  for  proofs.)  We  consider  only  this  case. 

Bifurcation  of  a steady  state  solution  of  (1.4)  can  only  occur  if  the  linear 

oiierator  1.  in  (1.5)  is  singular.  To  avoid  confusion  let  us  say  explicitly  that  L 
0 2 

operates  on  C ((0,t|,J(  ) with  homogeneous  Dirichlet  boundary  conditions.  Since  L 
2 

commutes  with  (S/Sf.)  , the  eigenfunctions  of  L may  be  sought  in  the  form 

fal 


(1.7) 


sin 


where  a,  b are  constants  and  2.  is  an  integer.  The  two  eigenvalues  of  L associated 
to  eigenfunctions  of  the  form  (1.7)  are  eigenvalues  of  the  matrix 


(1.8) 
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Zero  is  an  eigenvalue  of  (1.8)  if  and  only  if  B « where 
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(i.9) 
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Thus  the  first  bifurccation  jK^int  of 
(i.lO) 


+ 


(1.4)  is  given  by 

min  B,  . 
t 


For  most  values  of  the  {>arameters  in  the  problem  this  minimum  will  be  achieved  for 
exactly  one  value  of  and  then  the  first  bifurcation  will  be  from  a simpl.  igen- 

value.  However  if 

(1-11)  + 1)^ 

for  some  integer  k,  then  the  minimum  in  (1.10)  is  achieved  at  both  « k axvi 
* k + 1 , and  the  eigenvalue  is  doiible.  In  other  words,  when  the  zero  solution  of 
(1.4)  first  loses  stability,  it  loses  stability  simultaneously  with  respect  to 
disturbances  of  wave  nuinl>er  k and  of  wave  numf>er  k ♦ 1.  This  is  th«‘  case  w»-  study 
in  the  present  paper. 

Our  task  of  analyzing  the  bifurcationa  of  (1.2)  near  a double  eig**nval'je 

may  be  conveniently  divided  into  twt>  steps.  it  turns  out  that  (in  non-d«.*genei ate 
cases)  either  one  or  throe  now  solutions  bifurcate  from  the  trivial  solution  at  a 
double  eigenvalue,  and  various  comliinat ions  of  stability  or  instability  for  the  new 
solutions  are  fxissible.  In  all  there  are  five  different  cases.  Our  first  problem  is 
to  determine  which  case  occurs  as  a function  of  the  parameters  A,  D^,  D^-  (Only  two 
of  these  are  really  independent,  as  the  reriuirement  that  a double  eigenvalue  occur 
leads  to  the  condition  (1.11).)  The  second  step  inthe  analysis  is  to  study  the  effect 
of  various  perturbations  on  the  bifurcation.  There  are  two  natural  perturbations  to 
consider  in  this  connection.  One  is  to  change  A,  from  unperturbed  values 

which  satisfy  (1.11)  exactly  to  perturbed  values  which  satisfy  (1.11)  only  approximatelv . 
Observe  that,  no  matter  what  the  values  of  these  parameters,  (1.3)  still  provides  a 
spatially  and  temijorally  homogeneous  solution  of  (1.1),  (1.2).  The  effect  of  this 
f>er turbation  is  to  split  the  double  eigenvalue,  so  that  two  separate  bifurcations  from 
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the  trivial  solution  occur,  both  at  simple  eigenvalues.  This  perturbation  induces 


secondary  bifurcation  away  from  the  trivial  solution.  The  other  perturbation  we 
consider  is  to  replace  the  parameter  A in  (1.1)  by  the  function 

coshv'T  (C  - 

(1.12)  A(t)  = A„  , 

cosh/t  — 

wliicti  is  motivated  by  the  following  considerations.  In  the  derivation  of  (1.1)  it  is 
as.sumed  that  A measures  a chemical  concentration  which  is  fixed  by  the  experimenter. 

In  practice  concentrations  can  only  be  fixed  at  the  boundary  of  the  domain;  in  the 
interior  the  concentrations  will  be  determined  by  solving  a Ixaundary  problem 

, 2 

— - lA  = 0 on  (0,n) 

3C 

A(0)  = A(it)  = Ag  , 

which  lias  solution  (1.12).  Here  t measures  the  rate  at  which  A is  depleted  relative 
to  its  diffusivity.  Wlien  t “ 0,  (1.12)  reduces  to  the  previous  case,  but  when  r > 0 

an  explicit,  trivial  solution  of  (1.1)  analogous  to  (1.3)  is  no  longer  available. 

Indeed,  as  a result  of  this  perturbation  there  may  exist  a range  of  B where  (1.1)  has 
no  steady  state  solutions,  spatially  homogeneous  or  otherwise,  at  least  not  in  the 
neighborhood  of  (1.3).  The  effect  of  this  perturbation  is  similar  to  the  effect  of 
imperfections  on  bifurcation  at  a simple  eigenvalue.  (See  for  example  (2,9).) 

There  is  an  important  symmetry  present  in  (1.4)  that  restricts  considerably  the 
fxjssible  behavior  in  this  problem.  Namely  (1.4)  commutes  with  the  reflection 

(1.13)  Rw(e)  - w(ii  - C)  . 

In  other  words,  no  change  in  (1.4)  occurs  if  one  ma)(es  a coordinate  transformation 
t ' - n - t which  interchanges  right  and  left  endpoints  of  the  interval.  Observe  that 

R sin  a " (-1)  sin  It  . 

Thus  one  of  the  two  eigenfunctions  of  L at  the  double  eigenvalue  is  even  and  one  is 
odd.  The  two  perturbations  of  the  problem  mentioned  above  also  conmute  with  (1.13). 


f 


We  study  the  equilibrium  equation  associated  to  (1.4),  namely 

(1.14)  Lw  + N (w)  = 0 , 

by  means  of  the  Lyapunov-Scliroidt  reduction.  (See  §4  for  details.)  At  a double  eigen- 
value this  procedure  reduces  (1.14)  to  a system  of  two  equations  in  two  unknown.s, 
depending  on  a parameter.  Let  us  write  the  reduced  equations 

(1.15)  G(x,l)  = 0 
2 2 

where  G : F « F ]R  . Here  x = (x,y)  parametrizes  the  kernel  of  L and 
A * B - Bq,  where  B^  is  the  bifurcation  point.  Our  notation  in  (1.15)  does  not 
ifidicate  explicitly  the  other  parameters  in  the  problem.  We  will  show  below  that  at  a 

double  eigenvalue  the  reduced  equations  may  be  written 

2 2 
X + ny  + Ax 

(1.16)  G(x,A)  = 

cxy  + Ay 

where  c e F and  n = ±1.  (This  equation  holds  'o  cubic  terms,  and  in  non- 
degenerate cases  the  cubic  terms  may  be  transformed  -hange  of  coordinates.) 

The  symmetry  of  (1.14)  with  respect  to  the  transformation  is  reflected  in  the 

fact  that  the  first  component  of  (1.16)  is  even  with  respect  to  y while  the  second 
is  odd.  Thus  the  bifurcation  diagram 

(1.17)  {(x,A)  e F^  X » . g(x,A)  = 0} 
is  invariant  under  the  reflection  y ^ -y- 

The  first  step  in  our  analysis  is  to  compute  the  coefficients  c and  n in  (1.16) 
as  functions  of  the  various  parameters  in  the  problem.  It  is  then  a simple  matter  to 
determine  the  type  of  the  bifurcation  diaqram  by  invoking  the  classification  results 
of  (3).  (Strictly  speaking  to  obtain  the  stability  properties  of  the  bifurcating  solu- 
tions the  arguments  of  (3)  must  be  supplemented#  which  is  done  in  §§6  and  7 of  the 
present  paper.)  This  solves  our  first  problem#  that  of  describing  the  bifurcation 
diagrams  at  the  double  eigenvalue  proper.  For  the  perturbed  problems  we  again  refer 
to  (3)#  this  time  for  a result  that  an  arbitrary  small#  symmetry  preserving  perturbation 
of  (1.16)  may  be  described#  up  to  a certain  equivalence#  by  two  parameters.  In  more 


technical  ian^/uar^**,  an  of-r  i atf  universal  unfoHinq  of  (i.l^S)  r<*quirf*s  two  p'^rameters . 

(See  53  for  greater  detail.)  Th«*  fact  that  two  parameters  suffice  depends  strongly 
on  the  ; r^'sence  of  the  symmetry  (1.13)  - without  this,  five  parameters  would  be 
re<|uired  (21.  Tht;  two  parameters  of  the  universal  unfolding  may  be  identified  with  the 
two  perturbations  of  the  physical  problem  discussed  above.  The  fact  that  two  parameters 
suffice  for  th^  unfolding  means  Chat  additional  perturbat ions  of  (1.4)  - one  might 
for  example  treat  B as  a variable  concentration,  partially  depleted  in  the  interior, 
or  one  might  consider  boundary  conditions  other  than  (1.2)  - would  not  lead  to  more 
corolicated  behavior  than  that  already  obtainable  with  the  two  perturbations  considered. 

This  f)afjer  has  much  in  common  with  Keener  (Sal,  a reference  kindly  brought  to  our 
attention  by  W.  M.  Hay.  We  feel  that  our  paper  sheds  new  light  on  the  subject  on 
several  accounts,  '■(uite  apart  from  questions  of  rigor.  Perhaps  most  important,  the 
existence  of  a universal  unfolding  with  a known,  finite  number  of  unfolding  i>aram<>ters 
ensures  that  the  sort  of  parameter  exploration  undertaken  in  these  papers  is  a terminat- 
ing process  - although  new  complications  may  be  introduced  into  the  model  indefinitely, 
no  'lual itatively  new  behavior  will  result  after  the  parameters  of  the  universal  unfold- 
ing are  represented.  Secondly,  the  present  paper  emphasizes  the  imfxjrtance  of  symmetry 
in  tfiis  prol>l»'m,  without  which  a number  of  qualitatively  new  phenomena  could  occur. 

For  exampl*-  the  occurrence  of  secondary  bifurcation  depends  on  symmetry.  Consider  a 
j-rTturbat  ion  of  t>»e  r<*duced  equation  (1.16)  of  the  form 
(1.17a)  ' 

where  M IS  a 2 ' 2 matrix.  Such  pjerturbat ions  do  not  affect  the  existence  of  the 
trivial  solution  x ■ 0 but  (in  general)  do  split  the  double  eigenvalue.  It  can  be 
shown  that  if,  s<iy  c 0 and  n « tl,  then  secondary  bifurcation  occurs  for  the 
perturbed  {problem  if  ^nd  ^ M is  upper  triangular  (non-zero  diagonal  entries 

permitted).  Here  bifurcation  means  the  crossing  of  solution  branches,  not  merely  the 
existence  of  a limit  point.  The  first  perturbation  discussed  above,  namely  moving 
A,  Ij^,  slightly  away  from  values  satisfying  (1.11),  may  be  represented  in  the  form 

(1.17a),  where  however  the  fact  that  the  perturbation  preserves  symmetry  implies  that 


M IS  fiiaqonal.  In  particular  M is  upp<-r  trian'^ular , so  s«-.oiiju^/  ...i  ..  i.aLi.  Jo*'s 
o'.'cur.  Of  course  for  a symmetry  breakinq  perturbation  no  su''b  : ^rr:n  i.  r,  !>e 

♦'Xpected.  A final  difference  between  f Sa ) and  the  present  paper  is  t/).jr  w- 
the  two  distinct  perturbations  - Keener  does  not  considc^r  (1.12).  Prom  o\ir  j.  int  of 
view  it  is  riaturai  to  try  to  find  physical  representations  for  both  ot  the  unfoldinq 
parameters  in  the  universal  unfoidinq,  not  to  mention  the  considerable  physical  /.tercet 
of  (i.l2)  . 

The  following  notation  will  be  used  throughout  this  paper,  bet  0 - and 

”2 

let  A * 6 A.  Equation  (1.11)  may  be  re-written 
( 1 . IH)  A ^ k (k  + 1)D  , 

where  D without  a subscript  indicates  In  Figure  1.1  the  lines  labeled  fV,k  + 1) 

indicates  the  lines  in  the  (A,D)  plane  where  (1,18)  is  satisfied.  betwe«*n  the  lines 
(k  - l,k)  and  (k,k  1)  the  first  eigenvalue  is  simple  and  the  associated  r:>',  n- 
function  has  spatial  dependence  sin  kf.  If  (1.18)  is  satisfied,  the  first  oifuication 
occurs  when 

(l.ig)  B = (1  I Dp^) (1  + Dp^)  * 

2 2 
where  u^*=k  and  p2®(k  + l). 
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52.  Statement  of  results 


In  this  section  we  present  our  main  results,  deferrinq  the  proofs  for  later  sections. 
Illustrations  of  the  bifurcation  diaqrams,  as  defined  by  (1.17),  offer  tlie  most 
convenient  format  for  this  presentation.  Some  comments  about  Figures  2.2-4  below  may 
help  in  their  interpretation.  These  fiqures  are  intended  to  represent  three  dimensional 
bifurcation  diagrams,  the  three  coordinates  beinq  >,  the  bifurcation  parameter;  x, 
the  amplitude  of  the  eigenfunction  whose  profile  is  symmetric  with  respect  to  (1.13); 
and  y,  the  amplitude  of  the  anti -symmetric  eigenfunction.  The  orientation  of  these 
coordinates  shown  in  Figure  2.2  is  retained  throughout.  The  "Y"  or  diamond  shaped 
fig’jre  at  the  ends  of  some  of  these  bifurcation  diagrams  is  intended  as  an  aid  to 
visualization  and  does  not  represent  part  of  the  diagram  itself.  In  Figure  2.2  we  have 
shown  the  bifurcation  diagrams  of  the  unperturbed  problem  as  a union  of  straight  lines 
through  the  origin.  Here  *'T'*  denotes  the  trivial  solution,  which  coincides  with  the 
>-«ixis.  The  branches  labeled  T and  1 span  the  plane  of  symmetry  {y  = 0},  while 
branches  2 and  3,  when  present,  are  located  symmetrically  in  a plane  fcx  + > * 0} 
r^rpendicular  to  the  plane  of  symmetry.  (Cf  (1.16).)  In  the  perturbed  diagrams. 

Figures  2.3  and  2.4,  we  have  labeled  the  solution  branches  far  away  from  the  bifurcation 
{X>int  by  the  closest  solution  branch  of  the  unperturbed  problem.  In  all  cases  the 
perturbed  diagram  consists  of  two  conic  sections,  the  branches  T and  1 lying  in 
the  plane  fy  » O},  the  other  branches  lying  syirmetr ically  in  a plane  fcx  + 1 = const). 

In  our  bifurcation  diagrams  the  various  solution  branches  have  stability  assign- 
ments determined  by  the  spectrum  of  tiie  linearized  equation.  We  claim  that  at  most 
two  eigenvalues  of  the  linearized  equations  lie  in  the  unstable  (right)  half  plane  - 
this  follows  by  perturbation  theory  from  the  observation  that  the  unperturbed  problem 
has  a double  eigenvalue  at  zero  and  the  remainder  of  its  spectrum  lies  in  a half  plane 
fc  i ^ ) » where  f > 0.  In  the  fiqures  we  indicate  the  three  cases  of  zero, 

one,  or  two  eigenvalues  in  the  right  half  plane  by  labels  +s,  - , +U/  respectively. 

Here  the  sign  gives  the  Leray-Schander  degree  while  the  latter  distinguishes  between 
the  stable  and  unstable  cases  with  positive  degree. 


We  should  caution  the  reader  that  our  representation  of  the  bifurcation  diagrams 

is  purely  local.  It  is  known  fll  that  for  sufficiently  small  A the  solution  of 

(1.4)  is  unique  and  that  for  X bounded  the  solutions  satisfy  an  a priori  estimate. 

Thus  the  bifurcating  solutions  must  turn  around  in  the  large,  as  indicated  in  Figure  2.5 
for  the  unperturbed  case  Iq  (notation  defined  below). 

First  we  consider  the  unperturbed  case  when  (1.18)  is  satisfied  exactly.  As 
mentioned  above,  either  one  or  three  non-trivial  solutions  of  (1.4)  may  bifurcate  from 
the  trivial  solution  at  the  double  eigenvalue.  If  one  new  solution  bifurcates,  it 
may  be  either  unstable  or  stable,  while  if  three  such  bifurcate,  either  zero,  one,  or 
two  of  them  may  be  stable.  We  refer  to  these  cases  as  1^,  ^^^2 

respectively  - the  Roman  numeral  indicates  the  number  of  bifurcating  solutions  and 
the  subscript  the  number  of  them  which  are  stable.  It  turns  out  that  which  case  occurs 

depends  on  A and  D but  not  on  0 . The  type  of  the  bifurcation  as  a function  of 

these  two  parameters  may  be  determined  from  Figure  2.1,  where  we  have  identified  two 
sets  of  five  regions  in  the  A,  D plane  in  which  different  behavior  obtains.  The 
two  different  partitions  correspond  to  k odd  or  even.  Equations  for  the  boundaries 
of  these  regions  are  given  in  §4.  We  have  sketched  in  Figure  2.2  bifurcation  diagrams 
for  each  of  the  five  cases.  (In  this  paper  we  do  not  consider  the  degenerate  cases 
when  (A,D)  lies  on  the  boundary  between  two  regions  in  Figure  2.1.) 

We  now  consider  the  effects  of  changing  the  parameters  A and  D,  the  first  of 
our  two  perturbations.  (As  above  the  bifurcation  diagram  does  not  depend  on  0.)  If 
(A,D)  varies  along  one  of  the  lines  (k,k  + 1)  in  Figure  2.1,  no  qualitative  change 
in  the  bifurcation  diagram  will  occur  unless  (A,D)  crosses  into  a different  region. 


(This  is  of  course  a local  statement  - the  size  of  the  neighborhood  in  which  it  is  valid 
decreases  as  the  boundary  is  approached.)  Thus  only  one  of  these  paraimeters  has  a 
qualitative  effect  on  the  nature  of  the  bifurcation  diagram.  Let  us  take  6,  the 

change  in  D.  Making  6 non-zero  splits  the  double  eigenvalue  into  two  simple  eigen-  « 

} 

j values,  causing  secondary  bifurcation.  We  illustrate  the  effect  of  this  perturbation 

i 

; I 
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II 
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3 2 

Figure  2 .26;  Unperturbed  III^ 

for  the  five  cases  above  in  Figure  2.3.  The  diagram  depends  on  the  sign  of  6 and  the 
parity  of  k as  indicated.  Note  that  the  X-axis,  representing  the  trivial  solution, 
is  a part  of  all  these  diagrams.  Of  course  the  effect  of  the  perturbation  is  strongest 
near  the  center  of  the  diagram  - although  the  perturbation  is  uniformly  small  through- 
out the  interval,  the  unperturbed  problem  is  rather  singular  near  the  double  eigenvalue 
proper  and  hence  more  sensitive  to  perturbations.  Perhaps  the  most  noteworthy  feature  of 
these  diagrams  occurs  in  Figure  2.3bj.  (Let  us  remark  that  the  circle  in  this  diagram  lies 
in  a plane  perpendicular  to  the  two  lines.  Similarly  for  Figure  2.3a.)  The  solution 
branch  that  originates  from  the  first  bifurcation  exists  only  for  a small  interval 
above  the  bifurcation  point  before  it  is  reabsorbed  by  a secondary  bifurcation.  A 
similar  phenomenon  occurs  in  Figure  2.3e2.  It  is  also  worth  remarking  that  the  first 
bifurcation  can  be  either  super-,  trans-,  or  sub-critical. 

As  to  the  second  perturbation,  we  have  sketched  in  Figure  2.4  the  bifurcation 
diagrams  which  result  from  taking  f > 0 in  (1.12).  The  outcome  depends  on  the  region 
and  on  the  [)arity  of  k in  a somewhat  confusing  manner.  For  example,  the  unperturbed 
diagrams  when  (A,D)  belongs  to  region  1 or  5,  k odd  or  even,  are  all  the  same, 
namely  akxiive  if  6^0  the  perturbed  diagram  depends  on  the  parity  of  k but  not 

the  region,  while  here  if  E > 0 the  diagram  depends  on  the  region  but  not  the  parity  of  k. 
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values  of  f wouUi  lack  physical  si^^nif icanee.  As  it  happen.';,  however, 
for  the  various  unf>orturbed  diaqram.s  of  type  III,  the  diagram  which  wrjuid  result  if 
f *■  0 occur  anyway  for  (A,D)  in  a different  region  and  k of  the  opposite  parity. 

This  duplication  does  not  occur  for  the  cases  or  We  have  none  the  les.s 

included  the  diagrams  for  t 0 with  the  label  "unphysical",  since  these  diagrams 
could  well  appear  if  the  perturbation  were  caused  by  some  other  mechanism  without  a 
definite  sign,  as  for  example  imj-x^sing  boundary  conditions  slightly  different  from  (1.2). 
We  do  not  analyze  such  possibilities  hero,  but  only  mention  that  no  matter  how  complicated 
the  perturbing  mechanism  (assumed  symmetric),  the  diagram  which  results  is  determined 
guai  itatively  by  the  values  of  the  two  parameters  in  the  universal  unfolding  of  (l.lf.). 
(See  (2)  concerning  non-symmetr ic  perturbations.)  This  consideration  gives  the  diagrams 
here  a greater  significance  than  otherwise  apparent. 

By  far  the  most  interesting  of  these  diagrams  is  that  of  Figure  2.4e^.  (We  caution 
the  reader  not  to  try  to  imagine  a surface  spanned  by  solution  branches  I and  2 or  1 
and  3 - rather  think  of  2 and  3 as  lying  in  a plane  cx  + > * 0.)  Here  the  two  non- 
trivial solution  branches  of  positive  degree  change  from  unstable  to  stable  as  > is 
increased,  without  ever  encountering  a zero  eigenvalue  - in  other  words,  they  undergo 
a Hopf  bifurcation!  As  mentioned  above  a Hopf  bifurcation  can  sometimes  occur  from  the 
trivial  solution  in  ♦-he  unperturbed  problem,  but  only  at  much  larger  values  of  B; 
indeed  the  Hopf  bifurcation  of  Figure  2.4e^  can  occur  for  parameter  values  where  the 
unperturbed  problem  does  not  admit  any  Hopf  bifurcations.  Because  the  eigenvalues 

must  bfj  close  to  zero,  the  period  of  the  associated  limit  cycles  will  be  large, 

-1/2 

sf.»ecif ically  0{i.  ).  In  a subsequent  publication  we  will  consider  the  stability  and 

the  domain  of  existence  of  these  limit  cycles. 

Finally,  let  u.s  suppose  that  both  6 and  f are  non-zero.  Then  there  exist  six 
regions  of  the  4,  e plane  where  the  perturbed  diagrams  exhibit  different  structure, 

as  indicated  in  Figure  2.6  for  the  case  TTIq.  Note  that  regions  1,  2,  4,  and  5 in 

2 

the  figure  only  contain  points  <4,c)  for  which  c **  0(6  ).  We  feel  that  regions  3 and  6 
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Figure  2.5:  Behavior  in  the  large 

must  be  regarded  as  the  generic  cases,  since  points  ({,£)  whose  components  are  of  the 
seime  order  will  belong  to  these  regions.  The  diagrams  corresponding  to  regions  3 and 
6 are  those  illustrated  in  Figures  2.4Cj^  and  2.4c2<  respectively:  the  diagrams  correspons 
ing  to  the  four  thin  regions  are  required  to  effect  the  transition  between  regions  3 
and  6.  It  is  quite  possible  to  Imagine  situations  where  these  transition  diagrams  are 
relevant,  but  we  suspect  that  their  detailed  structure  is  too  specialized  for  the  model 
at  hand.  We  refer  the  interested  reader  to  (3]  for  a more  complete  discussion. 


E 


Figure  2.6 


§3.  Review  of  the  classification  results  of  [ ] 

In  the  next  section  we  will  apply  the  Lyapunov-Schmidt  reduction  to  the  equilibrium 
equation  (1.14)  and  will  obtain  a systfsn  of  the  form 


2 2 
a,x  + a.y  + 


b^lx  “ 0 


I a^xy  t b^ly  - 0 , 

at  least  modulo  cubic  terms.  These  equations  represent  the  most  general  reduced  equa- 
tions at  a double  eigenvalue  for  a system  governed  by  a quadratic  non-linearity  and 


[ossessing  the  synmetry 


where  J is  the  matrix 


JG(Jx,X)  - G(x,X)  , 


In  this  section  we  recall  some  of  the  definitions  and  results  of  (31  concerning  the 

classification  of  such  problems. 

2 2 

Let  G,H  : K » V D be  two  reduced  bifurcation  problems,  defined  near  the 

origin.  We  shall  call  G and  H contact  equivalent  if  there  exists  an  invertible, 

2 

2x2  matrix-valued  function  ^ and  a dif feomorphism  on  R x R of  the  form 

(x,X)  e*  (p  (x,X)  ,A  (A) ) such  that 

(3.3)  //(x,A)  = • G(p(x,A)  ,A(A))  . 

We  are  primarily  interested  in  bifurcation  diagrams  such  as  (1.17),  and  it  is  clear 
that  the  bifurcation  diagram  of  G is  not  changed  at  all  by  multiplication  by  an 
invertible  matrix.  The  dif feomorphism  (p,A)  represents  only  an  inessential  change 
of  coordinates  in  the  problem,  which  will  not  change  the  qualitative  nature  of  the 
bifurcation  diagram.  These  remar)(8  are  intended  to  motivate  our  definition.  (In 
general  transformations  such  as  (3.3)  can  change  the  stability  assignments  of  the  various 
solution  branches  in  a bifurcation  diagram.  However  the  degree  of  a solution  branch 
is  well  defined  since  we  require  that  and  dp,  the  differential  of  p,  have 

positive  determinants.  We  Ignore  the  stability  issue  for  the  time  being,  returning  to 
it  in  57.) 
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The  above  definition  neqlocts  the  symmetry  of  (3.1).  The  iat<*  coordinate 


l 


■I 


transformations  (0,A)  which  preserve  symmetry  must  satisfy  (3.2)  and  the  matrices 

T m»ist  satisfy 
x> 


(3.4) 


J ^ T 1 J 1 

JX,A  XX 


(In  practice!  (T.4)  simply  mi-ans  that  the  diagonal  entries  of  t are  even  functions 
of  y;  the  off  diagonal  entries,  odd.)  We  shall  call  two  problems  G and  M 
cquivariantly  contact  equivalent  if  (3.3)  is  satisfied  and  all  functions  have  the 
appropriate  symmetry.  Usually  we  abbreviate  this  phrase  to  equivalent. 

It  is  shown  in  I31  that  the  classification  of  problems  of  the  form  (3.1),  up  to 
e'^ui valence,  def^ends  only  on 


b a 

(3.5)  c « and  n “ sign(a,a„)  . 

Vi  ^ 2 


In  other  words  (3.1)  is  equivalent  to  (1.16),  where  c and  n are  given  by  (3.5). 

The  type  of  the  bifurcation  diagram  as  a function  of  the  parameters  c and  g is  given 
in  Figure  3.1.  Strictly  spcaiting  any  two  problems  of  the  form  (1.16)  with  different 


c 0 

0 < c < 1 

c > 1 

n * +1 

I, 

I„ 

HI 

— J 

0 

1 

n « -1 

1 

I 

0 

. _J 

2 

1 

Figure  3.1 


values  of  c are  inequivalent,  but  in  practice  the  qualitative  features  of  the  bifurca- 
tion diagram  are  unchanged  if  c remains  within  one  of  the  intervals  of  Figure  3.1. 
Moreover,  provided  c avoids  the  boundaries  of  these  intervals,  the  higher  order  terms 
tliat  were  neglected  in  writing  (3.1)  may  in  fact  be  transfontnd  away  by  an  auspicious 
choice  of  T , p , A . 

Reference  [3)  also  considers  perturbations  of  (3.1),  In  terms  of  t)ie  following  speci  f ic 

2 2 

definition.  By  an  unfolding  of  a bifurcation  problem  C s F » R « F we  mean  a 
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2 I 2 

smooth  m.i|.  f : K " t'  " such  that  F(x,A,0)  - G(x,>).  Thn  !.  variables  in 

the  third  argument  of  F will  he  called  unfoldirKj  parameters.  We  show  in  |3|  that 

'22 

X + ny  + lx  + e 

) (i.f))  F(x,>  ,u,f(,y)  = 

(c  + Y)xy  + (1  + 8)y 

is  a universal  unfoidinq  of  (1.10)  relative  to  equivalence,  assuming  c * 0,1.  This 
means  in  particular  that  given  an  unfolding  of  (1.16)  of  the  form  F(-,',t)  = G + tH, 
there  exist  smooth  functions  u(f).  Bit),  ylt),  defined  for  small  r,  such  that 
'^'0  ' ’ e<iuivalent  to  F ( - , ■ ,u  (t)  , P,  ( t ) , Y (t) ) . See  (31  for  a more  thorough  discussion 

of  ttiese  issues. 

Tlie  reader  will  note  that  our  universal  unfolding  (3.6)  contains  three  unfolding 
parameters,  while  we  indicated  above  that  two  parameters  suffice  to  describe  jierturba- 
tions  of  (1.16).  The  explanation  of  this  discrepancy  is  as  follows.  The  parameter  c 
determines  the  overall  character  of  the  bifurcation,  but,  as  noted  above,  once  it  is 
restricted  to  one  of  the  intervals  in  Figure  3.1,  small  changes  in  this  pareimeter  do 
not  affect  the  qualitative  nature  of  the  bifurcation.  Thus  we  regard  Y as  an 
inessential  parameter  for  our  present  purposes. 

The  following  table  indicates  the  correspondence  between  the  notation  of  the  present 
paper  and  that  of  (3)  in  naming  the  various  cases  which  arise. 


5_4 . The  I-yapunov-Schmifit  reduction  for  the  unperturbed  case 

Let  us  beqin  our  discussion  with  some  comments  about  the  interF^retation  of  (1.4) 

and  Its  associated  ecF\j i 1 ibr ium  equation  (1.14).  We  consider  (1.4)  as  an  evolution 

0 2 

er^ation  on  the  Banach  Sfiace  X ^ C ({0,ff),F  ).  We  could  of  course  regard  L as  an 
unbounde*d  linear  operator  on  X with  domain 

D = {wf  C^([0*tt1,1R^)  : w(0)  ^ w(Tr)  = 0}  , 

but  w(?  {>refer  to  regard  L as  a bounded  linear  operator  L : D •*  X » where  D has 

2 •* 
the  C topology.  Moreover  N,  considered  as  a map  from  D into  Xt  is  C 

Krechet  differentiable  - indeed  N is  already  smooth  considered  as  a map  from  X into 

00 

itself.  Thus  the  right  hand  side  of  (1.4),  namely  L + N,  defines  a C mapping 
t I D * X.  Note  that  which  is  a Fredholm  operator  of  index  zero. 

Lo*’  X^  = ker  L and  X^  ~ range  L.  We  shall  assume  below  that  X admits  the 
de'':om|>osi  tion 

(4.1)  X ^ ^ X^  . 

0 1 

This  means  that  all  the  generalized  eigenfunctions  of  L associated  to  the  eigenvalue 

2 

zero  already  belong  to  ker  Lj  in  other  words,  ker  L « ker  L,  We  also  suppose  that 
the  sF^ectrum  of  L restricted  to  is  properly  contained  in  the  stable  half  plane,  say 

(4.2)  (3(LlAr^)  c (j;  : RHi;  < -e} 

whern  t > 0.  Since  X C D,  it  follows  that  D ^ X ^ D , where  D ^ D X • Also, 

0 0 111 

L ■.  D ^ • X ^ is  a linear  isomorphism  between  the  two  Banach  spaces.  Let  E^,  be 
the  projections  associated  to  the  decomposition  (4.1). 

In  the  Lyapunov-Schmidt  reduction  one  eliminates  all  but  a finite  number  of 
components  of  w in  (1.14)  by  inverting  the  non-singular  part  of  L.  Specifically 
define  a mapping  W : 0^  implicitly  by 

(4.3)  E^'»(x  + W(x))  = 0 

for  X f X It  follows  from  the  implicit  function  theorem  that  (4.3)  is  soluble  in 
some  nelqhtxjrhood  of  zero.  Let  ’ ^0  ^ ^0  defined  by 

(4.4)  G(x)  - Eq*(x  + W(x))  . 
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Then  X ♦ W(x)  is  solution  of  (1.14)  if  and  or»Iy  if  G(x)  = 0,  ^ind  t'Vry 
of  (1.14)  has  this  form. 

For  'jomputations  it  is  convenient  to  have  a reformulation  usirri  oord  i r»a » er.  of 

the  l.yapunov-Schmidt  reduction,  carryinq  explicitly  ♦h*'  bifur'. jMon  pfir  im*  » * r /.  :.*  ♦ 

2 

be  a basis  for  X^.  If  we  param^ftrize  X^  by  P via  r he  rar  ✓ ► 

then  we  may  rewrite  (4.3)  as 

(4.5)  />W(x,>)  ♦ EjN(xw^  + yw^  ♦ W(x,>))  0 . 

* * * * 

bet  w^ , w^  be  a basis  for  ker  1.  , where  L is  the  adjoint  of  f.orrt-i,  /<  l . • i 

(4.6)  <w.,w.)*'i'^.. 

1 3 i) 

for  some  positive  constant  a.  (Here  and  below  < •,*)  denotes  th<;  l/'  . nner  produ'd.) 
Then,  in  t<  rms  of  the  coordinates  on  X^  above,  has  the  repre  .e^r -i f i or. 

* * 

(4.7)  E^w  » ( (w^ ,k >, (w^ , w>)  . 

2 2 

Thus  equation  (4,4),  defining  the  reduced  mapping  G : ' p » P , m »/  t.»-  r'  wr»**«  r« 

in  coirifjonents  as 

* 

(4.H)  G.(x,X)a(w.,'t>  (xw^  t yw-  + W (x , > ) ) ) , 

i - 1 1 2 - 

where  i » 1,2.  We  shall  always  take  w^^  to  be  even  with  rf*sf-<;ct  to  (I  .13)  -ukI 

* * 

odd;  w^  and  inherit  the  corresponding  parities  by  (4.6). 

We  may  perhaps  avoid  confusion  if  we  introduce  notation  to  indi'^^ate  th<'  d^'{'f*nd'  ncf* 
on  B that  is  implicit  in  L,  Let  us  write  Z<  « + AM,  where  is  giv'-n  by 

(1.5),  with  B fixed  at  the  bifurcation  point,  and 


Using  the  fact  that  » 0 we  may  rewrite  (4,8)  as 

* 

(4.10)  G.  (x, A)  - < w , (AM  + N)  (xw  + yw  + W(x , A ) ) > . 

—X  » 1 1 « •” 

We  now  start  the  computation  of  the  coefficients  in  (3.1)  as  functions  ot  t.)i»' 
various  parameters  in  our  problem,  Tn  the  following  lemma  d^N  denotes  th»*  second 
Frechet  derivative  of  N,  a bilinear  mapping  X ^ X X , 
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I.fiinina  'I . i : Th»*  fificofxi  or']t:r  'li-r  i vat  iv#;s  of  th<*  map  •'i  at  ori'Hn  a 


fjiverj  by 


* 1 * J 

- ( w N (w  . ,w  ) ) 
'Jx  . Jx  ^ j } 9 

J ». 


1 

- ■ < w ,Mw  , > . 

-lx  , 1 ' 3 

J 

Proof:  It  follows  by  d if f or i.-nt iat ion  of  (4.10)  that  (when  A = 0) 


(4.11)  . » < w,  .ri'^lKw,  t - , w + --  )>  + <w.,dN  ■ 7 > . 

-Ix^SXj,  1 3 -ix^  }.  -)Xj  1 

However  dK  is  zero  at  the  origin.  Thi.s  mean.s  that  the  second  term  in  (4.H)  vanishes 

f3W 

at  tht:  origin  and  that,  by  (4.5),  also  vanishes  there.  Tfiis  verifies  the  first 

j 

formula  in  the  l^mma,  and  the  argument  for  the  second  is  similar.  The  proof  is  complete. 

Wc!  do  not  reproduce  the  details  of  the  computation  of  the  coefficients  in  (3.1), 
but  only  record  the  landmarks  as  a guide  to  the  reader.  It  follows  from  (1.6)  that 
at  the  origin 


d^N(w  ,w  ) =»  2{^  u u + A(u  V + v u )] 
12  A12  12  12 


The  eigenfunctions  of  L and  L at  the  double  eigenvalue  are  given  by 


Wj  (? ) » sin  1?C 


w^((;)  . sin  kK 


w^Cf.)  = sin()c  + DC 


-1  - Dp., 


^1  " 9^2  ■ ^ 


is  a correction  factor  required  by  the  normalization  (4.6).  The  reader  will  note  that 
if  k Is  odd,  then  w^^  in  (4.13)  is  even  with  respect  to  the  symmetry  (1.13).  .Supixjse 
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that  the  two  wave  numbers  at  the  double  eigenvalue  are  t and  I + 1 . If  ?.  is  odd, 
we  shall  use  (4.13)  with  k - if  I is  even,  we  shall  use  (4.13)  with  k « -(f  * 1). 

By  this  ruse  we  can  always  arrange  that  the  smaller  wave  number  (in  the  algebraic- 
sense)  at  the  double  eigenvalue  is  odd.  In  this  way  we  may  observe  the  parity  conten- 
tion introduced  above  and  still  avoid  carrying  two  sets  of  formulas  for  the  eigenfunctions. 

On  substitution  of  (4.9),  (4.12),  and  (4.13)  into  the  formulas  of  Lemma  4.1,  one 
finds  after  appropriate  manipulation,  the  following  formulas  for  the  coefficients  in  (3.1). 

'2 


(4.15a) 


'1 


-t-2-^ 

li,  1 + Du, 


1 


-3  ■ 2'l 


(4.15b) 


(4.15c) 


(4.15d) 


(4.15e) 


Here  we  have  used  the  notation 


Du. 


2 “^2 

*2  ■ ^ I + DuY 


I, 


Du. 


Du, 


1 - 


1 + DU^  1 ♦ DU2 


b^  - f 


fi  2 «du2 


“ SAB  / sin^kf.  dC 


1 - BAB  / sin  kt  8in^(k  + DC  d5  . 

^ 0 

Note  that  I^,  I 2 have  the  same  sign  as  k. 

According  to  the  results  of  [3)  quoted  In  53  of  the  present  paper,  the  qualitative 
type  of  the  bifurcation  diagram  associated  to  (3.1)  can  change  only  If  one  of  the 
following  equations  Is  satisfied. 


(4.16) 


(a)  a. 


(c)  a. 


(b)  a. 


(d)  b^a3  - bjSj 


I>»t  us  consider  the  simplest  of  these  equations,  (4.16c).  On  multiplying  (4.15c)  by 
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(i  ♦ 'lividirK}  Ly  a r»on“Z#*ro  factor  w<‘  may  r^'Writ^*  M.l^yC)  as 

2 

i - f-*  ^^^2  “ ^ * 

-2  2 

W<‘  m.iy  comijino  the  definitions  of  §Z  to  show  that  A » D from  wliich  it  follows 

that  (4.if>c)  is  satisfied  if  and  only  if  A = 1.  The  reader  will  note  that  this  line 
aj  peiirs  as  one  of  the  dividing}  curves  in  Figure  2.1,  both  cases  (a)  and  (b)  . 

Consideration  of  equations  (4.iha)  and  is  similar.  First  invoke  (4.15) 

to  show  that  thf?sc'  equations  are  satisfied  if  and  only  if 


(a)  DU2  * i 


(b)  Du^  - 1 


respect! vely,  and  then  use  the  definitions  of  §1  to  express  (4.17)  in  terms  of  the 
parameters  A and  D.  If  k is  positive,  the  result  is  that  (4.if>a  or  b)  is  satisfied 


if  and  only  if 


(a)  D » (A  - 1)  , A < 1 


(b)  D =*  (A  - 1)  , A > 1 


resf/ectively;  if  k is  ntqative , the  formula  D = (A  - 1)  remains  valid  but  the 

inequalities  in  (4.1H)  are  reversed.  The  reader  will  again  note  that  the  parabola 
2 

0 » (A  - 1)  occurs  as  a dividing  curve  in  Figure  2.1.  The  remaining,  unidentified 
dividing  curve  in  the  figures  is  of  course  associated  with  (4.16d).  Its  equation  may 


iXi  written 


2f2(l  - D - (1  + Dp^)(l  - DU^)  = 0 , 


f .1/1  , __jJK  i)-_  , 

2 2'  1 (k  + 2)  (3k  + 2)  ' 


blit  (4.1'l)  doos  not  sefdn  to  admit  a simple  representation  as  a function  of  A and  D. 

In  Fiqure  4.1  we  have  tabulated  the  signs  of  the  coefficients  (4.15)  in  the  various 
regions  usln<j  the  notation  (3.5).  The  type  of  the  bifurcation  diagram  may  be  determined 
by  ':omparjson  with  Figure  i.l. 


V -ft;  V. 


55.  Lya^^unov-Schmidt  calculations  for  thf*  pfrrturhf-d  probl^^ms 

It  may  pjjrhaps  fjiv'-  a clearer  expxDsition  if  we  indicate  how  to  compute-  tl»e  unfold- 
inq  f/arameters  in  (3.6)  for  a general  perturt-at ion  of  the  equilibrium  equation  fl.l4)/ 
although  for  the  problciin  at  hand  we  are  able  to  avoid  the  most  unpleasant  aspf-cts  of 
this  calculation.  Supf>ose  (1.14)  is  modified  by  a perturbation  term  to  read 

(5.1)  i-Cw)  + F.p(w)  = 0 . 

l/*t  fc;^,  be  defined  as  in  §4  to  be  the  projections  relative  to  the  decomposition 

(4.1) ,  this  data  being  computed  only  for  t * 0.  We  define  the  reduced  equations, 

depending  on  c,  as  follows.  bet  W i F * be  defined  implicitly  by 

(5.2)  E^('t  tP)  (X  * W(x,(.) ) = 0 


(5.3)  G(x,£)  « EQ('t  + tP)  (X  + W(x,t))  . 

Translation  of  these  expressions  into  coordinate  notation  poses  no  special  difficulties, 
althouqh  it  should  be  mentioned  that,  unlike  in  §4,  the  function  W plays  a role  here. 
It  follows  from  the  natural  analoque  of  (4.5>  that  at  the  origin 

(5.4)  P , 

3e  1 

where  L ^ is  the  generalized  inverse  of  L.  Suppxise  we  define  coefficients  for  the 
l^erturbed  problem  by 

2 

G.(x,c)  -G.(x,0)  = cfa.  + y B..X.  + + hot  , 

1 1 1 13  3 lO 

3=1  •>  J 

2 2 2 

where  the  higher  order  terms  include  terms  of  order  f.  , r.x  , c>  . A straightforward 
calculation  shows  that 


a . * < w . ,P > 

X 1 

B..  =<w*,dP{w.)  - d^N(L'^E,P,wJ  > 

1]  1 ] 1 j 


«io  = If  - • 


The  only  significant  difference  between  these  formulas  and  those  of  bemroa  4.1  is  the 


appearance  of  a second  term  in  the  inner  products  of  (5.5b  and  c) , and  this  difference 


is  a direct  consequence  of  (5.4). 


Let  us  first  compute  the  unfolding  parameters  for  the  perturbation  associated 
to  (1.12).  It  is  readily  seen  that 

tA 

A(e)  = Aq  — ~ 5(1'  - £)  + . 

note  that  A appears  only  in  the  first  equation  of  (1.1),  and  there  only  additively. 
(Remar)c:  In  the  passage  to  (1-4)  we  still  linearize  about  a constant  function. 


namely  X = Y = B/A^.)  Thus  in  the  notation  of  (S.l)  we  should  define 


P = - J AQf.(Ti  - f)  i 


a i,erturbation  which  is  independent  of  w.  We  compute  from  (5.5a)  that  = 0,  as 
required  by  symmetry,  and  tliat 

2A 

(5.7)  (1  + DUj)  . 

)< 

The  computation  of  the  first  order  terms  would  be  rather  tedious  because  of  the 

presence  of  L ^ in  the  relevant  formulas.  However,  we  recall  from  53  that  the 
zeroth  order  coefficient  a in  (3.6)  dominates  the  first  order  coefficient  B,  provided 
they  are  of  the  same  order,  which  is  the  case  here,  since  (5.7)  shows  that  a is  non- 
vanishing. Thus  the  effect  of  the  perturbation  (5.6)  is  determined  by  the  sign  of  a 
(negative)  and  the  sign  of  the  coefficients  in  Figure  4.1.  The  reader  may  consult  13] 
to  checl<  the  validity  of  the  diagrams  of  Figure  2.4. 

The  perturbation  associated  to  changes  in  the  diffusion  coefficient  D admits  the 
representation 


1 0 ,2 

o k 


0 9 1 3? 


which  depends  linearly  on  w.  We  could  easily  substitute  (5.8)  into  (5.5)  and  evaluate 
the  resulting  expressions  - this  perturbation  differs  from  the  preceding  one  in  that 
(lere  P vanishes  at  the  origin,  eliminating  the  troublesome  terms  with  L However, 

there  is  a direct  way  to  assess  the  effect  of  (5.8),  which  moreover  sheds  some  insight 
on  the  problem.  This  perturJoation  splits  apart  the  double  eigenvalue  but  does  not 
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affect  the  existence  of  the  trivial  solution.  We  r;.  see  fron  (?.^))  that  there  is 


essentially  only  one  way  to  do  this,  namely  to  chan‘:<e  6 but  to  keep  a equal  to  zero. 
Thus  to  describe  the  perturbed  diaqram  we  need  only  know  which  mode  bifurcates  from 
the  trivial  solution  first,  and  we  can  decide  this  by  an  independent  argument. 

Let  us  recall  (1.9),  the  fo*7»ula  foi  the  bifurcation  point  of  the  mode  with  wave 
number  i from  the  tiivial  solution 

2 

(5.9)  B = 1 + — - + De  + r . 

eor 

2 2 

Observe  that  we  may  write  (5.9)  as  a function  of  c real  variable  Dfi.  , say  * f (Di  ). 
By  hypothesis 

tnin{f  (D£^)  : I = 1,2,3,...} 

1 

is  assumed  at  two  distinct  integers  k and  k + 1#  as  indicated  in  Figure  5.1.  It 
may  be  seen  by  inspection  that  increasing  D lowers  B^^  and  raises  other 

words,  increasing  D makes  the  mode  with  smaller  wave  number  (in  absolute  value) 
bifurcate  first.  This  is  the  behavior  portrayed  in  Figure  2.2.  To  facilitate  the 
reader’s  checking  this  statement,  we  mention  that  the  solution  points  of  the  mode 
with  odd  wave  number  are  located  in  the  plane  of  symmetry  of  the  bifurcation  diagram, 
while  those  with  even  wave  number  occur  in  symmetric  pairs. 


/ 


! 


1 

I 


§6.  On  the  relation  between  the  Lyapunov-Schmidt  and  center  manifold  reductions 

The  Lyapunov-Schmidt  reduction  enumerates  all  rest  ix>ints  of  (1.4)  but  provides 
no  information  about  the  dynamical  behavior  of  this  equation.  In  this  respect  reduc- 
tion of  (1.4)  via  the  center  manifold  theorem  is  preferable,  although  somewhat  less 
straightforward.  We  begin  this  section  with  a brief  review  of  the  projJerties  of  the 
center  manifold,  referring  to  (4)  or  (81  for  details  not  provided  below.  We  continue 
to  use  the  notation  introduced  in  54  for  the  Lyajjunov-Schmidt  reduction. 

The  center  manifold  M is  a finite  dimensional  submanifold  of  X parametrized 
by  X e X^.  More  precisely 

(6.1)  M = [x  + V(x)  : X e Kp  close  to  zero) 

where  V : K -*  D , satisfies 
0 1 

(6.2)  (E^  - dV  - E^ltlx  + V(x))  = 0 . 

The  reduced  mapping  H : -►  in  the  center  manifold  context  is  defined  by 

(6.3)  ll(x)  ■=  Eq4(x  + V(x))  . 

Thus  (x  + V(x)  : H(x)  = 0}  provides  an  alternative  enumeration  of  the  solutions  of 
(1.14).  At  the  same  time,  however,  the  trajectories  of  (1.4)  tend  to  trajectories  of 
the  ordinary  differential  equation 

(6.4)  ||  = H(x) 

in  the  following  sense.  Let  U be  an  appropriately  small  neigliborhood  of  zero  in  X, 
let  w(t)  be  a solution  of  (1.4)  such  that  w(t)  e U for  0 t ^ T,  and  let 

«(t)  - ||EjW(t)  -V(EpW(t))l|  , 

the  distance  between  w(t)  and  its  projection  into  W.  Then 

4(t)  < (0) 

for  0 t ^ T,  where  t is  defined  by  (4.2)  . The  constant  C depends  only  on  U in 
the  above  data.  Unfortunately  this  result  does  not  state  in  general  to  what  extent 
the  projection  x(t)  EjjW(t)  is  approximated  by  a solution  of  (6.4),  but  if  4(0)  -0 
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(i.e.,  if  w(0)  e M)  , then  x(t)  is  actually  a solution  of  (6.4).  Moreover  if  S C 
is  an  attracting  set  for  (6.4),  then  (x  + V(x)  : x e S)  is  attracting  for  (1.4). 


Equation  (6.2)  expresses  the  condition  that  the  flow  direction  iKw)  be  tangent 
to  M when  v e M.  To  see  this  let  us  define  a smooth  mapping  f X -<■  by 

f(w)  = Ej^w  - V(EqW)  , 

so  that  W = f ^(0).  We  may  write  the  condition  of  tangency  as 


^ f (w(t) ) =0  . 
dw 

since  M is  contained  in  D,  — is  well  defined  for  v e M.  Application  of  the  chain 

at 

rule  leads  to  (6.2) . 

In  one  sense  (6.2)  is  a small  perturbation  of  (4.3);  namely,  since  dV,  ♦,  and 
Eq  • d*  all  vanish  at  the  origin  the  perturbing  term  may  be  expected  to  be  small  near 
zero.  On  the  other  hand  the  new  term  in  (6.2)  involves  derivatives  of  V,  which  malces 
an  existence  proof  for  (6.2)  by  a direct  perturbation  argument  problematic.  Indeed 
this  existence  question  is  discussed  in  [41  in  terms  of  the  non-linear  sanigroup  of 
transformations  generated  by  ♦.  Moreover  because  of  the  singular  nature  of  this 

go 

perturbation,  (6.2)  does  not  necessarily  admit  C solutions.  Although  there  are  c 
solutions  for  arbitrarily  large  )t,  the  size  of  the  domain  of  existence  decreases  as  )t 
increases.  However,  the  fact  that  dV,  ♦,  and  • d*  all  vanish  at  the  origin  does 
have  the  following  consequence:  in  computing  derivatives  of  V at  the  origin  from 
(6.2),  the  second  term  will  always  contribute  to  lower  order  than  the  first.  Thus 
for  example 


(6.5) 


,-1.2.. 

- -L  d ♦(e.,e^) 


where  L : K -*■  0 is  a pseudo-inverse  and  (e.)  is  a basis  for  the  same 

11  1 0 
2 

formula  holds  for  3 W/3x.3x.. 

1 3 

There  is  a minor  technical  issue  to  be  addressed  before  the  results  of  (4)  are 
formally  applicable.  We  must  modify  (1.4)  outside  of  a neighborhood  of  zero  in  order 
to  verify  the  hypotheses  of  [4] . Let  x be  a C function  with  compact  support  on 
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th»*  litiit**  M jm»*ns lorj.i I .i«  r ^ ^ nc*ar  z»‘ro.  W**  nKxiify  (1.4)  to  rf'ad 

=*  /.(w)  * x(K^w)N(w)  , 

SC)  tli.it  tli<-  o<pi,it  ion  IK  lin»‘-ii  whnn  K^w  is  larqe.  It  is  then  a siin[)le  matter  to 
v*T»fy  the  e‘.t  im.it«*s  n**efl»Ti  m (41  to  construct  the  center  manifold.  For  this  vrjrifica- 
t ion  the  re.ider  should  not#*  the  following  fact:  N(w)  vanishes  to  second  order,  so 
hy  restrictino  rlw  stjf/f  irt  of  x we  may  arrange  that  xN  and  its  first  derivatives 
arc*  .IS  small  is  may  be  necessary.  Also  if  depends  smoothly  on  a finite  numbfjr  of 

; ar.un*  t"r‘; , if  is  shown  in  (4|  fhat  the  center  manifold  depends  C -smoothly  on  those 
{<aram‘  t er  s . 

Th«-  fact  that  the  center  manifold  is  only  finitely  differentiable  would  lead  to 
certain  technical  compl icat ions  if  we  attf^pted  to  apply  the  center  manifold  theorem 
{lirectly.  The  followind  lemma  provides  a relation  between  the  center  manifold  and  the 
I.yapunuv-Sclim idt  reductions  that  may  be  exploited  to  avoid  these  complications.  The 
matrix  in  the  lemma  depends  on  the  parameters  B,  A,  in  the  problem  but  wo 

do  not  indicate  this  dependence  explicitly. 

h.i;  There  exists  a C matrix-valued  function  o such  that  the  reduced 

- X 

mappings  0 and  !1  of  formulas  (4.4)  and  (6.3)  are  related  by 

(h,h)  M(x)  - o • G(x) . 

X 

2 

Moreov» r n is  invertible  and  n * I + Oix  ). 

X X 

iT^^of : As  a technical  device  we  introduce  additional  parameters  into  the  arguments 

of  G and  II.  bet  us  define  an  extended  function  ‘P  : D ^ X by 

•Kwjy)  a ♦(w)  + y . 

W«>  r»d#fine  all  the  functions  above  to  include  these  extra  parameters,  indicating  the 
change  with  a twiddle.  Thus  for  example  ^ defined  by 

(6.7)  (E^  - d^V  • Eq)*(x  + V(x,y))  - - y - 0 , 

wh'Tf  have  used  the  relation  E^y  * 0 to  discard  a term.  It  turns  out  that 
(6.H)  G(x,y)  » G(x)  + y . 
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It  is  clear  from  this  equation  that 


(6.9)  f (x,y)  : G(x,y)  = Ol 

is  a smooth  submanifold  of  * X^,  which  is  the  reason  for  introducing  the  additional 
parameters . 

The  fundamental  observation  in  this  proof  is  that  G(x,y)  =0  if  and  only  if 
= 0 - both  reductions  yield  all  rest  points  of  the  equation.  Therefore  each 
component  of  //(x,y)  vanishes  on  the  smooth  manifold  (6.9),  so  by  Taylor's  theorem 
each  such  component  may  be  written  as  a linear  combination  (with  variable  coefficients) 
of  the  functions  defining  (6.9),  namely  the  comjDonents  of  G itself.  In  other  words 
we  have 

(6.10)  W(x,y)  = o G(x,y) 

xy 

for  some  matrix-valued  function  o.  We  obtain  (6.6)  on  restricting  (6.10)  to  y « 0. 

The  proof  that  = I O(x^)  is  based  on  the  conputation  of  various  derivatives  of 

(6.10)  indicated  below.  Given  this,  it  follows  that  a is  invertible  for  small  x. 

X 

The  derivatives  we  shall  need  first  are 


t 


(6.11) 


3G 

9x 


0,  d G = I,  = 0 

1 y 3x.3y^ 


for  X = y = 0.  The  last  two  relations  are  obvious,  in  view  of  (6.8),  and  the  first 
follows  from  (4.4)  on  observing  that 


(6.12) 


Eq  • d*  - 0 


We  claim  that  H satisfies  the  same  relations  (6.11).  Now 
(6.13) 


//(x,y)  = E^'Kx  + V(x,y))  + y . 


The  first  two  relations  in  (6.11)  follow  inmediately  from  (6.13)  by  differentiation, 
if  (6.12)  is  recalled.  As  to  the  third  relation  in  (6.11),  we  have 


(6.14) 


3V 


3x^3yj 


again  ma)cing  use  of  (6.12)  to  drop  a term.  We  argue  that  3V/3y^  = 0 as  follows. 


Differentiate  (6.7)  with  respect  to  and  evaluate  at  x = y = 0.  This  yields  simply 
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57.  Stability  of  the  bifurcating  solutions 

Thn  goal  of  tho  present  section  is  to  verify  the  stability  assignments  made  in 
drawing  the  bifurcation  diagrams  of  52.  We  do  this  by  ad  hoc  arguments  which  def.end 
strongly  on  the  symmetry  (1.13). 

The  bifurcation  diagram 


(7.1)  { (x,A)  ! H(x,l)  ■»  0} 

enumerates  the  equilibrium  solutions  of  (1.4).  Here  H is  the  reduced  mapping  on  the 
center  manifold,  as  defined  by  (ft. 3).  In  the  discussion  below  we  suppress  the 
dependence  of  H on  all  external  parameters,  including  the  bifurcation  parameter  1. 
According  to  56  the  stability  or  instability  of  any  equilibrium  solution  of  (1.4) 
enumerated  by  (7.1)  is  determined  by  the  stability  or  instability  of  the  corresjmnding 
rest  fioint  of  the  ordinary  differential  equation  (6.4).  However  by  l.emma  6.1  we  may 
express  H in  terms  of  G,  the  reduced  mapping  of  the  Lyapunov-Schmidt  reduction, 
and  G in  turn  may  be  expressed  in  terms  of  the  universal  unfolding  (3.6), 


(7.2) 


G ■ T F » p . 
x 


On  combining  these  observations  we  see  that  (6.4)  may  be  written 

— “ O T F • p (x)  . 
dt  XX 

If  we  Introduce  the  change  of  coordinates  x*  • p(x),  we  may  compute  that 

(7.3) 


^ - v(*')  • 


where  i ■ dp  • o • t . To  summarize,  we  saw  in  53  that  the  rest  points  of  (1.4) 

XXX 

may  be  enumerated  by  the  zeros  of  F;  here  we  see  that  the  stability  properties  of 
these  rest  points  may  also  be  obtained  frcm  F through  analysis  of  (7.3). 

The  stability  or  instability  of  rest  points  of  (7.3)  is  of  course  determined  by 
the  real  parts  of  the  eigenvalues  of  the  Jacobian  of  this  equation,  which  at  a test 
point  of  (7.3)  equals  stability  assigisnents  on  the  bifurcation  diagrams 

in  52  were  made  by  an  inspection  of  the  eigenvalues  of  dF,  or  in  other  words,  were 
made  according  to  the  stability  of 

5^-F(x)  . 


1 


(7.4) 
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J 


Thus  our  t<»sk  in  pr^'sont  section  is  to  show  that  the  real  [>arts  of  the  eigenvalues 
of  dK  are  not  changed  by  mu  1 1 ipl ica tion  by 

The  basic  fact  Wf*  will  use  in  this  derivation  is  that  ^ satisfies  (3.4),  which 

X 

foll<iWs  from  the  fact  that  each  of  the  factors  in  the  definition  of  t satisfies  (3.4). 

X 

Thus  jn  the  plane  of  symmetry  (i.e.,  when  y = 0) , is  diagonal.  Wo  claim  in  fact 

that  the  diagonal  entries  of  ^ are  both  [jositive  in  this  plane.  ft  suffices  to 
look  at  the  origin,  since  is  non-singular  and  cannot  have  a vanishing  diagonal 

entry.  At  the  origin  = I and  may  be  ignored.  Taking  an  appropriate  mixed  deriva- 
tive* of  {1.2)  we  find  that  at  the  origin 

• -li)  • ' 

only  one  non-zero  term  results  from  differentiation  on  the  right  since  F vanishes  to 


a#'corid  order.  however  it  may  be  seen  from  (3.6)  that  ^ ^ 


diagonal  matrix  witli  entries  b^,  b^  given  by  (4.15),  both  positive.  Thus  the  diagonal 
entries  of  the  product  * dp  are  [Xjsitive  at  the  origin,  and  being  diagonal,  these 

matrices  commute.  This  proves  the  above  claim.  In  the  following  we  shall  omit  the 
bar  in  i , as  we  have  no  further  occasion  to  refer  to  the  individual  factors  in  this 

X 


If  tfje  parameters  c and  n in  (3.6)  have  the  same  sign,  we  claim  that  (7.3) 
and  (7.4)  admit  a common  Pyapunov  function,  and  hence  have  the  Scime  stability  properties. 
The  Pyapunov  function  is 

2 , x^  n(l  +3)  2 

(7.5)  't(x,y}  * — ♦ nxy  + > -:p  *■ S Y ♦ ux  , 

i 2 c 

where  c*  * c t y.  We  suppose  that  Y is  small  so  that  c*  and  n also  have  the 
riamr-  sign.  To  check  that  (7,5)  is  a Lyapunov  function  for  (7.4),  we  compute  that 

(7.^.)  <F,  grad  ♦ > * (x^  + ny^  + Ax  + (i)^  + ^ (c'xy  + (A  . B)y)^  . 


( F , grad  ^ r.  | P | 

where  f • mind,  2n/c*),  so  ^ does  indeed  increase  along  the  orbits  of  (7.4).  For 
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(7.3)  we  must  compute  <T  F,  qrad  (P  > . Let  us  write  r “ + t',  where  is  the 

X X 0 X 0 

vaiue  of  at  the  origin.  Now  is  a jxasitive  definite,  diagonal  matrix,  and 

It  is  easily  seen  that 

(7.7)  ^ '7’^  ~ ' 

for  some  appropriately  small  positive  fi . On  the  other  hand  ^ > vanishes 
to  higher  order  at  the  origin  and  may  be  dominated  by  (7.7)  in  a sufficiently  small 
ne ighlaorhood . This  proves  the  above  claim  and  thereby  verifies  the  stability  assign- 
ments made  for  the  three  cases  1^,  7IIj-  (Incidentally  since  (7.4)  admits  a 

Lyapunov  function,  no  Hopf  bifurcation  is  possible  for  these  cases.  Cf.  below.) 

It  remains  to  verify  the  stability  assignments  for  the  two  cases  and  IIl2> 

in  which  c > 0,  n •"  0.  (The  case  c 0,  n > 0 does  not  occur  for  the  problem  at 
hand.)  The  bifurcation  diagram 

f (x,l)  : F(x,l)  - 0) 

consists  of  two  conic  sections,  one  in  the  plane  of  symmetry  y » 0 and  one  in  the 
plane 

(7.8)  c'x  + 1 + B ” 0 . 

We  refer  to  these  as  the  symmetric  and  asymmetric  solutions  respectively.  The  stability 
properties  of  the  symmetric  solutions  present  no  problem,  because  in  the  plane  of 
symmetry  and  dF  are  diagonal  matrices,  the  entries  of  being  positive;  thus 

the  eigenvalues  of  dF  and  of  • dF  are  real  and  have  the  same  signs.  Our  analysis 

of  the  asymmetric  solutions  is  based  on  the  following  two  facts,  proved  below.  In 
these  statements  large  only  means  close  to  the  boundary  of  an  appropriate  neighborhood 
of  zero. 

(7.9)  In  case  III2,  the  asymmetric  solutions  are  unstable  for  sufficiently  large 
negative  1 and  stable  for  sufficiently  large  positive  X. 

(7.10)  At  most  one  Hopf  bifurcation  is  encountered  along  the  asymmetric  branches 
an  A is  increased  from  large  negative  values  to  large  positive  values. 


I 
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AssumirMj  (7.’i)  and  (7.10)  thfi  readcir  may  verify  our  assi<^nm<jnts  of  stability  in 
the  remaining  dia^jrams  of  52  by  tracinq  the  various  solution  branches  in  from  iarqe 
values  of  A,  usinq  the  principle  of  exchanqe  of  stability  at  bifurcation  points.  Th^- 
essential  idea  here  is  that  the  stability  assiqrunent  alonq  a smooth  branch  of  the 
bifurcation  diaqram  can  only  chanqe  at  a Mopf  bifurcation  i>oint,  and  (7.10)  provides 
the  necessary  control  to  show  that  there  is  a unique  assignment  ;x3ssibie.  For  example 
in  Figure  2.46^  no  llopf  bifurcation  is  |xjssible,  while  in  Figure  2.4e^  a Hopf  bifurca- 
tion is  retniired.  It  should  be  mentioned  that  Figure  2.4bj^  was  drawn  assuming 
1 c < 2;  a slightly  different  diagram  results  if  c > 2»  and  we  exclude  the  degenerate 
case  c ~ 2 from  consideration  below. 

Both  (7.9)  and  (7.10)  follow  from  the  same  estimates.  On  using  (7.8)  to  eliminate 
X from  tlie  first  e(|uation  in  (3.0)  we  find  that  the  asymmetric  solutions  lie  on  the 
conic  section 


(7.U) 


I J_ 
,2  ” c* 


Now  along  the  asyiwnetric  solutions 


dF 


2x  t X 


-2y 


(c*y  0 J 

2 

Thus  dftt  dF  » 2c*y  ^ 0,  so  both  eigenvalues  have  real  parts  of  the  same  sign;  the 
same  statement  holds  for  • dF,  and  we  may  determine  the  sign  by  examination  of 
the  trace.  A trivial  calculation  yields 

(7.12)  tr  T • dF  ■ T (2x  t A)  t T,„(cy)  + t.,  (-2y)  . 

X 11  12  21 


We  remifKl  the  reader  that  and  are  odd  functions  of  y,  so  that  the  last 

2 

two  terms  in  (7.12)  are  0(y  ).  If  we  use  (7,8)  and  (7.11)  to  eliminate  x and  y 
from  (7.12)  we  find  that 


(7.13)  tr  T^dF  - (1  - ~)X  + hot) 

2 

whure  terras  that  are  of  order  \ , a,  or  B are  called  higher  order,  as  the  first 
term  In  (7.13)  will  dominate  the  others  if  the  parameters  are  appropriately  restricted. 
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In  case  ^^^2'  have  0 c 1,  so  that  the  coefficient  of  ^ in  (7.13)  * 

is  negative.  Thus  tr  t • dF  is  positive  for  large  negative  indicating  eigen- 

i 

values  with  positive  real  part  or  instability  as  claimed  in  (7.9);  and  similarly  for 

large  positive  >.  A Hopf  bifurcation  can  occur  only  if  (7.13)  vanishes,  which  will 

* * 

happen  for  exactly  one  value  of  >,  say  A . If  1 yields  real  solutions  y in 
(7.11)  the  associated  bifurcation  diagram  will  possess  a Hopf  bifurcation  point;  other- 
wise not. 
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